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^ ; ABSTRACT 

The chiral expansion of the low energy processes n° — > 77 and r\ — > 77 is reconsidered with 
particular emphasis on the question of the evaluation of the two low-energy parameters from 
£(6) Z wri ich are involved at chiral order six. It is shown how extensive use of sum-rules and 
saturation with resonances as well as constraints from asymptotic QCD effectively determine 
their values. Predictions for the widths are presented for both standard and non-standard 
values of the quark mass ratio m s /rh. A precise relation is established between the usual 
phenomenological n — n' mixing parameters and those of the chiral expansion. The large size of 
the chiral correction to the r? decay can be understood on the basis of a simple counting rule: 
0(1/ N c ) ~ 0(m q ). It is shown how this counting rule eventually allows one to include the n' 
into the effective lagrangian in a consistent and systematic way. 
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1. INTRODUCTION 

The 77 electromagnetic decay is a neat low-energy process, free of any strong final-state 
interaction, so one would expect that the chiral expansion should converge as fast as it seems 
to do in the case of the masses. Yet, at leading order, the width is predicted to be too small by 
a factor of three compared to experiment. The problem of the 0(p 6 ) corrections to r\ — > 27 (and 
7T° — > 27) has been addressed several times in the literature|T|[0](who computed the chiral loop 
contribution) ||[|J. Strictly speaking, chiral perturbation theory (xpt) is unable to make any 
prediction in this case because there are precisely two low-energy (LE) parameters, which are 
finite and called t\ and t[ below, which cannot be determined at present from other low-energy 
data (in principle the amplitudes for 77 — > 7r°7r°77 or 77 — > 37r° could be used as they contain 
just the same parameters). However, shortly after the xpt has started to be developed in a 
systematic way || || along with its list of LE parameters, a great deal of progress was made 
in the art of relating, in a rather precise way, the size of these parameters to the properties of 
the low-lying resonances 0@. In particular, as || have emphasized, whenever one can write 
down a rapidly converging dispersion relation it is natural to expect that saturation from the 
first low-lying resonances should provide an adequate approximation. 

In this work, we intend to exploit in this manner specific convergence properties of dispersion 
relations, which can be shown to hold in QCD in the chiral limit, in order to estimate the two 
LE parameters which occur in the electromagnetic decay of the it and the 77. These sum- 
rules are similar to the ones that were established for Lio0 an d for L7|10| on the basis of 



two-point functions. The parameters £1 and t[ are related to three-point functions instead. We 
will find out that t[ is proportional to the square root of — L 7 and to the amplitude A(r]' — > 
77). This amplitude must be taken from experiment and this information feeds back into 
the LE parameter, as is usual for the properties of the light resonances 0. Similarly, the 
parameter ti may be shown to encode information on the 7r(1300) resonance electromagnetic 
decay|3[]. Fortunately, (because hardly anything is experimentally known on this decay channel) 
we show that this contribution is numerically dominated by a contribution from the asymptotic 
behaviour of the three-point function which turns out to be canonical. A similar behaviour was 
exploited previously in the case of three-point functions involving one scalar current and was 
claimed to provide an explanation, based on chiral symmetry, for the somewhat unexpected 
electromagnetic properties of the scalar mesons []TT . 



Besides the purely phenomenological application, once it is ensured that the LE parameters 
are estimated to a reasonable level of accuracy, one can address the question of the rate of 
convergence of the expansion in powers of the quark masses. The conventional determination 
of the light quark running masses m u and m^ [0 was recently challenged [pfl, based on an 



analysis of the violation of the Goldberger-Treiman relation. This analysis suggests that the 
value of the quark mass ratio 2m s /{m u + m ( i) is two to three times smaller than in the standard 
Xpt, i.e. r ~ ?"2 = 25.9. If true, this necessitates to rearrange the quark mass terms in the 



chiral perturbation series in a different way [14] [p~5f| |16[ (the so-called generalized xpt) which 



one would expect to converge more rapidly. In practice, distinguishing between the xpt and its 
generalized variant is, curiously, not so easy. It is likely that only very precise measurements 
of 7r — 7r scattering lengths could settle the issue of which one is correct. It turns out that the 
amplitude A(rj — > 77) is very sensitive to the value of r when r is in the vicinity of the value r 2 . 
Nevertheless, we will show that it is also possible to reproduce the experimental results under 
the assumption of a much smaller value of r. 

The plan of the paper is as follows. The next section contains the derivation of the sum- 
rules for the parameters t\ and t[, the main results are contained in the formulae (|1^) and (p3|) . 



Application to the 77 decay amplitude is then discussed in sec. 3. In particular, we establish a 
connection between the chiral expansion description and the phenomenological representation 
in terms of an 77 — 77' mixing angle which is widely employed (e.g. ||17|| ). The chiral corrections 



to the 7T° decay are also worked out here. This analysis suggests that the value of F^o quoted 
in the literature is incorrect. The description of the 77 decay amplitude in the generalized xpt is 
presented at the end of sec. 3. The 77 — 77' system is of interest also in connection with the large 
N c expansion since, in the large N c limit, the 77' is a ninth Goldstone boson. This suggests 
still another expansion scheme where 1/N C is considered as an expansion parameter together 
with the quark masses and the momenta. A natural choice is to count one power of 1/N C on 
the same footing as one power of a quark mass. This point of view is discussed in sec. 4 (in 
particular, one observes the appearance of two different mixing angles but only one of them is 
relevant to the 77 decay) and the results are compared to those of the standard chiral expansion. 

2. LOW-ENERGY PARAMETERS FROM A CHIRAL SUM-RULE METHOD 

The current-algebra prediction for the amplitudes describing pseudo-scalar meson decays 
into two photons is contained into the following term of the canonical Wess-Zumino lagran- 
gian|T§[g 

r WZ _ N c nvap , m 

(4) ~~ 32tT 2 F (PV^Vap > UJ 

where the totally antisymmetric tensor is normalized such that e 0123 = 1 and <f) = J2 <fii\ j i=l,8. 
The subscript in eq.(|lj) is a reminder of the chiral order. Part of the higher chiral corrections 
are contained in chiral lagrangian terms which are also proportional to the e tensor. At order 
six there are three independent terms which are relevant for our purposes (using relatively 
standard notations, see e.g. ref.0) 

£%? = ^ va "{ *i < x { - ] && > +t[ < x H >< /£>/# > 

-tt 2 <d x f^{f^\u u } >+...} (2) 



As discussed in [fZC], all the other potentially relevant terms that one can write down can be 
reduced to the above three by the use of the so called Shouten identity and the equation of mo- 
tion. Among the constants appearing in (|2|) only t 2 has been estimated previously pl||22|[23 . 



The parameter t 2 is phenomenologically interesting in that it controls the corrections to decays 
like 7T° — > 77* which were discussed in f22| . This contribution vanishes when both photons are 



on-shell. In this case, the terms t\ and t' x are the only contributions from CW, Z to the processes 
7T° — > 27 and 1] — > 27. These two parameters are finite, reflecting the fact that the loop contri- 
bution is finite and, as a matter of fact, vanishes 0j@. The parameter t[ is analogous to some 
extent to the parameter L7 of the standard 0(p A ) lagrangian|6|j [|T0|j . It picks up a contribution 
from the pole of the rj' and, from the point of view of the large N c counting, it is of order 
0(N 2 ) instead of being 0(N C ) at most, like the other parameters. From this point of view one 
expects it to play a dominant role. The second parameter, t\ is unrelated to the 7/ resonance. 
In order to obtain an estimate for t±, t[ and £2 let us take the chiral limit and consider the 
vector-vector-pseudoscalar correlation function, i.e. 

d%w3PVrW(p 2 ,g 2 ,r 2 ) = Jd i xd 4 ye^ + ^ < 0\T(%(x)j b Mj e P (0J)\0 > (3) 

with r = — (p + q), a, b, c = 1 to 8 and where 

\a \c 

jfa) = ^(aOyT/^tX), Jp{x) = iip(x)—-f 5 ip(x), (4) 

are the vector and the pseudo-scalar currents respectively. In the case of the singlet pseudo- 
scalar current, i.e. when the index c is set to c = in (|3]) we define a function Hyy P from 
exactly the same formula. T^yvp differs from Hvvp because of the presence of the U(l) axial 
anomaly. 

Firstly, let us perform the low-energy expansion of Tlyyp. At order 0(p e ) one has to include 
a tree-level contribution from £%? ', a one-loop contribution involving one vertex from £|^ z and, 
finally, a tree contribution involving the parameters t\, t[ and £2 from C^i? (§). In the minimal 
subtraction scheme the scale dependence which arises from the loop contribution is cancelled by 
that of the parameter t 2 [fl)[22fl) the first two parameters are scale independent. The expansion 
is as follows 

IWp(p 2 , q\ r 2 ) = 2B J ^^ + 16^ + 4t 2 V -^- j + (chiral loop) (5) 

The first term in ([|) is the pole contribution from the canonical Wess-Zumino lagrangian. In 
the pseudo-scalar singlet channel, now, we no longer have a Goldstone boson pole contribution 
so that the chiral expansion starts with a constant term: 

U VVP {p 2 ,q 2 ,r 2 ) = 2B {m 1 + 48t[} + O(p 2 ,q 2 ,r 2 ) (6) 

For moderate values of the momenta p, q ~ 1 GeV it is natural to assume that Hyyp and 
Tlyyp are dominated by the low-lying vector and pseudo-scalar resonance poles. A very useful 
constraint arises in the chiral limit from the behaviour in the asymptotic regime p, q — > 00. 
Indeed, a simple calculation in QCD shows that one has0: 



2 The author is indebted to Marc Knecht for noticing a sign error in the first version of the manuscript 



In other terms, under a scaling of the momenta p —>■ Xp, q — > Xq, Hvvp scales as 1/A 4 if we 
let A go to infinity. The leading term in the asymptotic behaviour has the property of being 
canonical, i.e. it does not contain powers of logarithms caused by anomalous dimensions. This 
is because the scalar condensate < iptp > has the same anomalous dimension as the scalar (or 
the pseudoscalar) current and the vector current carries no anomalous dimension. In deriving 
(|7|) we have performed the operator product expansion a la SVZ|24|], i.e. taking into account 



the non-perturbative feature that the quark and gluon condensates are non-vanishing. SVZ 
argue that doing this enlarges the domain of applicability of the asymptotic behaviour down to 
the 1-2 GeV domain. It is therefore natural to demand that the two domains match smoothly. 
This is really useful only in the chiral limit because the correlation function does not blow 
up asymptotically and one need not resort to tricks like the Borel transformation which are 
necessary in the more general situation. Resonance saturation in the low to medium energy 
region amounts to approximate the correlation function by a meromorphic function having 
simple poles at the location of the resonance masses. Reducing to a common denominator, one 
can in principle have an arbitrary polynomial in the numerator. It is here that the asymptotic 
conditions come into play and limit the degree of the polynomial. In the present case, where 
the asymptotic behaviour is canonical, the term of highest order in the polynomial gets exactly 
determined. Let us for instance include one nonet of vector resonances in the vector channel and 
only the pion octet in the pseudoscalar channel. Then, we obtain a very simple representation: 

n ,222^ - B o F o (P 2 + q 2 + r 2 ) + a 

UvMp > q > r) = 2(p>-M v ^-M v y* (8) 

where My is the vector meson mass in the chiral limit. This construction bypasses the use 
of effective lagrangians for resonances and yields automatically the most general form of the 
amplitude containing the right poles and obeying the appropriate asymptotic constraints. In 
the case above one has a single arbitrary parameter, a. In writing down (||) we have ignored 
the contribution of the 7r(1300) (tt 1 ) multiplet. If we include the tt' pole contribution into the 
amplitude, we obtain a more complicated representation, 

, , , (~B F 2 (p 2 + q 2 + r 2 ) + a)(r 2 - M P 2 ) + b(p 2 + q 2 ) + cr 2 
nyyp(P ' g ' r) " 2(p> - M v 2 )(q* - M v *)r*(r* - M P *) () 

which contains two additional parameters b and c. Mp is the tt' nonet mass in the chiral limit. 
This representation reduces to the preceeding one in the limit where the mass Mp is sent to 
infinity. It is useful to rewrite (|9|) by separating out the various poles in r 2 

1 




2(p 2 - M v 2 )(q 2 - M } 



2^ 



X 



q 2 )(B F 2 + b/M P 2 ) c + b/M P 2 (p 2 + q 2 ) 1 

+ is — rrs B o F o ) (10) 



2 r 2 -Mf 



The parameters a, b and c may be related to properties of the resonances which will feed 
back into the low energy parameters which interest us upon expanding the correlation function 



around p 2 = q 2 = r 2 = and comparing with the representation (|5|). Here, we must exercise 
some care with respect to the chiral loop part which is not explicitly displayed in (H). Setting 
p 2 = q 2 = makes it vanish[l[]||. This allows us to unambiguously identify the parameter a 
from the residue of the pion pole: 

a = ^B M v i (11) 

and to find a relation between t\ and c: 

— 1 c 

32£ ii = z(B Fn + o) (12) 

2M/ V ° ° M P 2) V ; 

In the case of t%, we must consider the terms proportional to p 2 . The chiral loop generates a 
contribution proportional to p 2 /r 2 (which cancels out the scale dependence arising from £2) and 
furthermore generates nonanalytic terms which have no counterpart in the parametrization (^). 
We will content ourselves with a simple estimate of £2 valid in the leading N c approximation: 
in this limit, the chiral loop needs not be taken into account as it is subleading in N c . Ignoring 
the loop, one finds 

4B t 2 = — Lj(-% - B F 2 ? ) (13) 

There remains to find estimates for b and c. As usual, physical amplitudes are extracted 
from the Green's functions by taking the residue of the appropriate resonance poles and dividing 
out by the meson-current coupling constants. Using that, we can express c, to begin with, as 
a ratio of amplitudes: 

c = atane^f- -Hf (14) 

A(tt° ^77) 

where G is an angle which parametrizes the strength of the coupling of the tc' to the pseudo- 
scalar current 

< 0\] P (0)W b >= 5 ab B F t&nG (15) 



in the chiral limit (the rationale for introducing an angle here is explained in ref . [|TT|] . ) As 
will be recalled in the next section (see fl62f )) tanO can be related via a sum-rule to the LE 
parameter L 8 . The amplitudes A(P — * 77) where P is any pseudo-scalar meson are normalized 
everywhere in the text such that the width: 

T = ^M P A 2 (16) 

Similarly, there is a simple relation between the parameter b and the ratio of amplitudes A(ir — > 

pi)/M* ^77): 

A{-k ^ P1 ) f -2eF v \ M v 2 mi7^ b ^ C\7\ 

-1 + x, x = —{B F +—^) (17) 



A{tx -»• 77) V M v J a v u M f 

where Fy is the coupling of the vector field to the vector current. In eg. (|I7"D , x measures the 
deviation from the vector meson dominance principle (VMD). Exact VMD means that x — 0. 



Using the experimental value f\ T = 68 ± 7 KeV for the decay width p + — > 7r + 7 (together with 
F v = 150 MeV and M v = 770 MeV) we obtain x = 0.022 ± 0.051 which, as has long been 
well known, is fairly close to exact VMD. If we let Mp go to infinity in (|H]) we get a deviation 
from VMD of the order of 20%. For a better accuracy, we have to use a finite n' mass. In this 
description, VMD is realized from a cancellation between the parameter b and the asymptotic 
term BqFq (see ( |TT| ) ). Furthermore, assuming the cancellation to be exact ensures that the 
form factor associated with the matrix element < 7|j^|vr > satisfies an unsubtracted dispersion 
relation. 

We can now express the parameters t± and £2 in terms of experimentally accessible resonance 
properties: 

-1 



U 



F S + Nc(Mv)\ me A(* v 7 ) 



M v 2 4tt 2 \Mp) A(vr^77) 



UM V 2 

N 

i! = s# |1+I » (18 > 

Since the value of x is compatible with zero, the expression above for £2 reproduces the one 
obtained previously ||21|| p2| . Concerning t 1; the contribution proportional to A(tt' — > 27) is 



identical to the one which was identified in ref. || but the first contribution (which comes from 
the asymptotic term in the three-point function) was not included in that paper. We claim, on 
the contrary, that this contribution is likely to dominate. Experimentally, there exists an upper 
bound for the decay 7/(1295) — > 27 rate, T < 0.3 KeV. Assuming ideal mixing in the 7r(1300) 
multiplet together with the estimate tan0 ~ 2 (which is obtained from the sum-rule quoted 
in sec. 4 (^) which relates it to the value of the LE parameter L 8 ) this bound implies that the 
first term in t\ is larger than the second by a factor of at least three. This means that keeping 
the first contribution in ti, i.e. ignoring the n' contribution, should provide, at least, the right 
sign and the right order of magnitude for this parameter. 

Let us now consider the flavour singlet pseudo-scalar channel. We can write a representation 
of the three-point function analogous to (p~0|) : 

1 
2(p 2 - My 2 )^ 2 - M v 
[ a' - (p 2 + g 2 )(J? F 2 + b'/M P ,') d + &'/M P , V + g 2 ) 1 

1 ^Mj + ^M^ BoF °j (19) 

Here, Mp> is the mass of the singlet state in the 7r(1300) nonet in the chiral limit. At leading 
N c order we expect nonet symmetry to hold and thus Mpi = Mp. Furthermore, in this limit, 
the residues of the poles in ( |T9"D should be equal to the ones in flT0|), i.e a = a', b = b' c = c'. 
One can express t[ as a difference between Hvvp an d ^-vvp m ^ ne limit p 2 = q 2 = r 2 = 
(again the loop contribution vanishes): 

9QB t[ = lim \U° VVP (0, 0, r 2 ) - IL/yp(0, 0, r 2 ) + §^} (20) 



n u - : x 



3 All the experimental numbers are taken from the 1992 edition of the particle data book |2£ 



which gives 

96B £ = --L (-?-» + -^ - -%) (21) 

01 2M/W/ M F 2 M P 2 ' V ; 

Here, one expects a strong cancellation between the c and the c' term. In fact, the first term 
in the parenthesis is of order 0(N 2 ) while the sum of the last two terms is of order N®. It is 
easy to check that adding further resonances would also generate corrections which are of order 
N® and which are suppressed by inverse powers of the resonance masses. Therefore, it appears 
that retaining only the first term in ([H]) should constitute a rather solid approximation. Next, 
using ( |I9| ) we can express the parameter a' in terms of the amplitude A(rj' — > 27). Introducing 
the coupling of the rj to the pseudo-scalar singlet current 

<0\j° P (0)\ V '>=B o G v , (22) 

we obtain the following representation for the parameter t[. 

^ = 2^^^ (23) 

(the amplitude is again normalized as in ( |I6|) ) At this point we have all the necessary ingredients 
to discuss the chiral corrections to the decays 77 — *■ 77 and 7r° — > 77. 

3. APPLICATION TO tt° AND rj DECAYS 

Let us begin with the 77. Collecting the chiral lagrangian pieces up to order 0(p 6 ) one 

obtains: 

e'F^F^fo 1 5-2T2 1 

£ *- 2 ^ 1 6^F„^ l 1 + ^ r ' + (1 " r2)ri / (M) 

where we have introduced for convenience the dimensionless quantities proportional to t\ and 

t'v 

256vr 2 _ , . 1024vr 2 „ , , , 

Ti = — g— M^ r; = -^—M^ (25) 

In ( p4|) r2 is the quark mass ratio m s /m which must be expressed at 0(p 2 ) precision, i.e. 

r 2 = 2M 2 K /Ml - 1 ~ 25.9 (26) 

and, as usual F^ u = l/2e flual3 F al 3. In order to compute the 77 decay amplitude from (|24[) we 
must use that 0s = 0r? at chiral order 0{p 2 ) while 8 = 0^F o /F^ at order 0(p 4 ), ignoring the 
it — rj mixing which effect can be shown to be negligible for the 77 decay. 

Since the loop contributions vanish for on-shell photons U the decay rate is computed 
by simply using (^4j) at tree level. Using relation fl23|) for t[ one can write the amplitude for 77 
decay in the following way 

A(, - rr) = -A- {- + ^^1 + ^(r 2 - i)^4^t(„< - 77 ) (27) 



In this formula, we have replaced F everywhere by F n : this does not modify the first term, 
and in the others, this replacement amounts to a modification of order 0(p s ). One observes 
that the rf amplitude has crept in via the sum-rule formula for the LE parameter t' x . 

Using this form of the amplitude we can make contact with the traditional parametrization 
of the 7] and rf decay amplitudes in terms of an 77 — rf mixing angle. This will allow us to 
establish definite relations between the parameters appearing in this representation and those 
of the chiral expansion. The representation which is usually employed contains three parameters 
6» , A and A 8 (see e.g. ]T7||) 

. , . a / cos 6*o v8 sin 9 

A( v _► 77 ) = —=—— - 



V37iF n K A 8 A 







A(V' - 77) = * (^ + ^^o) (28) 



y/3nF w ^ A 8 A 

Let us eliminate Ao from the second equation and replace in the first one. We obtain a form 
exactly similar to (|27|) , which allows us to express A 8 and 9$ in terms of LE parameters: 

1 F - + 5 -^n Tl (29) 



cos 6»o A 8 F v 

The usual assumption made in the literature to identify A 8 with the ratio F n / 'F n is seen to be 
justified only to the extent that the correction from t± is negligible. Here, we find that this 
correction is of the order of 10%. Concerning the mixing angle, we obtain 

-^'--^f (30) 

A relation between the mixing angle and the parameter L 7 was first derived in || based on an 

analysis of the pseudoscalar mass formulas. We obtain a similar, but slightly different, formula 

upon using the sum-rule relation between Gy and L 7 (see (^) below). The two formulas do 

coincide at first order in the quark masses, as they should. 

Let us now make some numerical estimates. At order 0(p A ), first, one has the well-known 

result 

(4) _ a 2 M v 3 

1 ++*> ~ 192tt3 F, 2 " ( j 



(In numerical applications we use the value F n = 92.4 ±0.2 MeV||26||). At order 0(p ) now, we 
can use F v = (1.3 ±0.05)^ provided by and t\ is given by ([18]) and the following discussion. 
In ([27]) we still need to evaluate the coupling constant Gy. An estimate may be obtained from 
the sum-rule for the LE parameter LylTOf 



48M/ 



(32) 



This result is derived in the chiral limit under the only assumption that the 7/ pole is the 
dominant contribution in the dispersion relation. Using the numerical value of Lj, L7 = — (0.4± 

9 



0.15) 10~ 3 H one deduces that G v > = 133 ± 25 MeV up to a sign. In the large N c limit 
we expect to have Gy = F n which incites us to adopt the positive sign. Using this value 
in the expression (^) for the amplitude, together with the experimental value of the width 
T(j]' — > 77) = 4.29 ± 0.19 KeV, we obtain the rj width in the following form 

T(r] -> 2 7 ) = 172 [(0.77 ± 0.03) + 0.09 + (0.72 ± 0.15)1 2 = 430 ± 98 eV (33) 



where we have displayed the contributions in the same order as they appear in (p7|). One 
observes that the third contribution, which comes from the parameter t[ largely dominates over 
the one from ti, in agreement with the expectation from large N c arguments. The contribution 
from ti is not completely negligible, it increases the result by about 50 eV. If we had chosen 
the opposite sign for Gy (implying a huge deviation from the leading N c estimate) then the 
width would have been of the order of 3 eV. The central value is slightly below the experimental 
result (the average of the two-photon data at present give T = 510 ± 26 eV ) but the relatively 
large error bar on L7 generates an uncertainty of nearly 100 eV on our prediction (as a matter 
of fact, within this uncertainty, the result is compatible with both the old Primakov result 
(r = 324 ±46 eV) and the photon-photon results). The uncertainty on L 7 is tightly correlated 
with the uncertainty on the value of the quark mass ratio r, which is r = 25.7 ± 2.3 in the 
standard xpt at 0(p 4 )|J]. Roughly speaking, the upper bound on T corresponds to the lower 
bound on r and vice versa. 

Using our formula ( |30| ) with Gy derived from L 7 , we obtain for the mixing angle 

6 = -(18.4 ±3.6)° (34) 

which is smaller than the one quoted in [ID], #0 = — (22±4)°. We should note here that the value 



of Lj that we used does not take into account the recent work which have provided estimates 
for the corrections to Dashen's theorem [ 27 1. Part of the large error bar in L-j is generated from 



the uncertainty associated with these corrections. 

Let us now turn to the tt° decay. In contrast to the rj case it is crucial here to take the quark 
mass difference m^ — m u into account. The lagrangian including the 0(p 6 ) corrections reads 

e 2 F lw F^ v cb'3 ( , 5 rrifj — m n . m m„ — rrid m ,1 , N 

^3-77 = ^ 2F 1 + (1 " o ■ T i + U — d T[ 35 

1671^ F I 3 m d + m u m u + m d J 

Clearly, one expects the chiral corrections here to be much smaller than in the 77 case because 
they do not involve the strange quark mass. The wave-function renormalization, in this case 
can be written as follows at 0(p 4 )|| neglecting terms which are quadratic in e 



At 0{p 2 ) one has to set e\ = e% 




(36) 



V3m d -m u 
4 m s — m 

10 



(37) 



and, numerically, e ~ 1.00 10~ 2 and e 2 — 1.11 10 -2 . Expressing next t[ in terms of the mixing 
angle 9 and the rf decay amplitude as before, we can write the amplitude for the pion decay 
in the following form: 



A(w -> 77 ) = — — 1 + -?= + (1 + r^-em - 3e tan# A(r}' -> 77) (38) 




One notices that F^o appears here and not F n = F n + which is experimentally well determined. 
Theoretically, there are two contributions to the difference between F n o and F n +. One arises 
from the quark mass difference nid — m u : a simple estimate (neglecting chiral logarithms) gives 

F*° = F^(l - |(^ - l)e 2 + \t\ - ei e 2 ) (39) 

which shows that this contribution is negligibly small, of the order of 10~ 4 . The other contribu- 
tion is purely electromagnetic and, being of order 0{clqed) — 10~ 2 it should be the dominant 
one. Unfortunately, no estimate of this effect seems to be available at present. 

Using the value of 9 as discussed above, one can express the rate including the 0{p & ) 
corrections as 



-r = 7.73 



1 + (6.41 - 2.49 + 12.5)1(T J = 7.98 ± 0.08 eV (40) 



In the parenthesis we have displayed the contributions from the terms proportional to 62 , T\ 
and tan 9q respectively. The value of the error includes the uncertainty on F n as well as that 
from L-j. Comparison with the experimental result, V = 7.74 ±0.60 eV, suggests that the value 
of F n o should be somewhat larger than F n . Despite the rather large error bar on the ir° width, 



the value quoted in ref.||25||, F n o = 84 ± 3 MeV, is incompatible with our chiral expansion of 
the amplitude. 

Let us finally discuss the description of the amplitudes rj(n°) — > 77 in the generalized xpt 
approach. The 0(p 4 ) lagrangian in the Wess-Zumino sector is obviously the same as in the 
conventional xpt, being independent of the quark masses. The leading corrections are given 
by exactly the same two terms, proportional to t\ as t[ as before, which are counted as 0(p 5 ) 
rather than 0(p s ). The main difference then, is that the value of the quark mass ratio r is no 
longer determined from the 0(p 2 ) xpt but is left as a free parameter. As a function of this 
parameter, we can express the r\ amplitude as follows: 

a ( F 5 — 2r 1 

A( V -> 77) = -^r [Y n + ^- Tl ^| - tan ^o(r)A(r/ - 77) (41) 

Concerning the it , we note first that the mixing parameters depend on r, for instance 
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where A GMO = 3M V 2 /M n 2 — 2r 2 — 1 and e has the same meaning as before (|37|). The amplitude 
for 7T° decay can be written as 



10(1 -r) 5-2r 
[2e + e(r)] tan o (r) A(»/ -> 77) (43) 



^ 77 ) =_|i + _u + [i + _L r t e+ e(r )lr l( r 



The expression for e 2 (r) is rather lengthy and contains several 0(p 3 ) low-energy parameters 
which are not precisely known, however, an approximate expression can be derived 



e 2 (r ) = e(r) +e 



2, n ^r-r 2 ,F i 2 



P^i)^ff(^-i) («> 



The determination of the mixing angle angle 0q(t) goes in the same way as before. Instead 
of eq.flSOD, we obtain here 



n ,, V2, s 2mB Gr,> 1 M n / 2(r 2 - r)(r - 1) \ 2 , . 

tan W = -(r - 1)-^ = ^ (-A GMO + L2_JL. ) j (45) 



The second equality is obtained upon using the sum- rule (p2|). Instead of L7, one must use the 
corresponding LE parameter Zq in the generalized xpt, which appears at 0{p 2 ) and is given 
at this order from the deviation to the Gell-Mann-Okubo mass formula (see ||16|| ). In order 
to determine Ti(r) one starts, as before, from the sum- rule ( |l8l) for t\. The coupling tan0 is 
estimated by using, again, the sum- rule ( |62"D replacing L 8 by the corresponding parameter, A , 
in the Gxpt. One ends up with the following expression: 

-4rr 2 M n 2 F n 2 n M w ( XM S 2 - (1 - \fM^ 2 \* A(n> ^ 77) 2(r 2 - r) 

llrj " 3 M/ [ } M P { M P 2 -M S 2 ) A(7r^ 77 )' r 2 - 1 

(46) 

Let us now consider some numerical results in this formalism for a small value of r, r = 10 for 

example. Using the experimental bound on ^(1295) — > 77 we find 

T 1 (r = 10) = -6 10' 3 (0.68 + X), \X\ < 1.38 (47) 

where X is the term proportional to A(ir' — > 77) in the formula for T\(r). This term generates 
an uncertainty of the order of 30 eV in the value of the rj width and of the order of 0.1 eV in 
the 7T° width. The remaining contributions to the 77 width are shown below: 

r(?? -* 77) |r=io = 172 [0.77 + 0.02 + 0.95] 2 = 521 eV (48) 

where we have used that F^/F^ is practically the same as in the standard xpt and that the 
value of the mixing angle for r = 10, as obtained from (f45"|) , is #0 — 23.8°. This result is at the 
upper limit of the allowed range in the standard xpt. The various factors in (EST) are, however, 



widely different from the corresponding ones in (30): r << r2, IviiBq << M w , but this is 
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compensated by the third factor, Gy ~ 8F n which is much larger than in the standard xpt 
(and deviates more from the large N c result). The ir° width, finally, can be expressed as 



7T 



r(7T° -> 77)|r=io = 8.07 -^eV (49) 



Both the results for the ir° and the 77 fall in the range of values which are also allowed in the 
standard xpt. 

4. MIXED LARGE iV c AND CHIRAL EXPANSION 

The proper way to introduce the pseudo-scalar singlet field 4>q(x) into the chiral lagrangian 
was first discussed in PS |. The most general lagrangian at order 0(p 2 ) is discussed in ||. We will 



follow their method here, which consists in introducing a source 8(x) associated with the winding 
number density uj = (lG7r 2 )~ 1 TYG llu G flu . Under an axial U(l) transformation, exp(i/3(x)), one 
lets the source 8(x) transform to 8(x) — 2tr/3(x). Then, (apart from the canonical WZ part) the 
effective action is invariant. In the singlet sector, three kinds of invariant building blocks are 
available -D M 0o = d^o ~ 2tra M , D^O = 8^9 + 2tra M and the combination -\/60o/-^o + 9. For the 
purpose of discussing the spectrum and electromagnetic properties it is sufficient to consider 
quadratic polynomials in <po{x). The canonical WZ action is minimally extended using nonet 
symmetry. 

At order two the mixing angle is predicted to be 9p ~ — 10°[l|]f| which is too small by a 



factor of two compared to phenomenological determinations |T7j • It has long been recognized 
that going to 0{p A ) seems to improve the situation]]]]. In ref. |H| it is proposed to use the 0{p A ) 
lagrangian with the rf integrated out in conjunction with the 0(p 2 ) lagrangian which includes 
the rf. Although this procedure gives the correct magnitude for 9p one may wish to include 
the rf at 0(p 4 ) as well. Furthermore, the lagrangian proposed in contains terms which 
are leading in N c as well as terms which are subleading. Outside of the large N c expansion, 
including the rf explicitly into the chiral lagrangian does not make much sense since the rf mass 
is rather large, larger than the mass of the p or u mesons, for instance. On the contrary, this 
looks reasonable if one thinks in terms of the large N c expansion, since M 2 , is 0(1/N C ). This 
suggests to generalize the chiral counting in order to include the large iV c expansion: a natural 
choice is to count a factor which is 0(1/N C ) on the same footing as a quark mass factor[29~|. We 



will call this counting rule the mixed chiral expansion. It should be clear that this expansion 
is different from the usual chiral expansion. For instance, at leading order, the eta mass does 
not satisfy the Gell-Mann-Okubo relation and the ratio F n /F v differs from one. However, it is 
also a systematic expansion which could, in principle, be carried out at arbitrary high order. 
We first generalize the chiral matrix U to include a nonet of fields: 

C7 = [/expi(^M) (50) 

v Fq ' 

At order two of the mixed expansion one has the two terms which are 0(p 2 ) x 0(N C ) as well 
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as a single term which is 0(p°) x 0((N C )°): 

£ (2 ) = X {< D ^D»U + X ] U + W X > ~l M o(^y o + ^) 2 } (51) 

It is not difficult to extend this to the next order. We note first that chiral loops need not be 
considered since they are of chiral order four and they are suppressed by one power of N c , so 
we effectively count this kind of contribution as order six. Furthermore, the two terms of the 
standard 0(p 4 ) lagrangian: 

£ = L 6 < X j U + W X > 2 +L 7 < X f U - W X > 2 (52) 

are both suppressed by the Zweig rule here so they are also effectively of order six (the L's pick 
up contributions from the resonances which are not included in the chiral lagrangian, this is 
why L7 is different from L7 ). The only terms that we need consider are 

£ (4) = U < D^D»U( x m + W X ) > +L 8 < X ] U X ] U + fP X &X > (53) 

In addition to these terms, which are 0{p i ) x 0(N C ), we have to consider the terms which are 
0(p 2 ) x 0((N C )°). We can form two terms of this kind: 

F 2 
6 V '~F 

Furthermore, in the Wess-Zumino sector, we can form a term which is 0(p 4 ) x 0((N C )°) and is 
thus on a similar footing with the true 0(p 6 ) corrections: 

c ,wz = h ^ a , K f M f W > (V g^ + 0) (55) 

Before proceeding, we note that the LE parameters L 5 , L 8 and the quark mass ratio r = 
m s /rh have to be re-determined using the logic of the mixed expansion scheme adopted here. 
In particular, L5 and Lg are scale independent at order four now. Outside of the rj — rf sector 
one may simply use the formulae of and drop the loop contributions. L 5 , to begin with, is 
related to the ratio F K /F n ~ 1.22: 



£' ( 4) = kxDpfoDPfa + ik 2 -±(Vs^ + 0) < X ] U - W X > (54) 



F — F 



2 Z?2 

, - 2.31 1(T 3 (56) 

a 8(Mi-M|) l ; 

This value is very close to the result of L^(fi = M v ) = 2.31CT 3 which is not very surprising, 
since the kaon and the eta chiral logarithms are very small for this value of the scale and the 
pion chiral logarithm, being proportional to the pion mass squared, is small anyway. 
Next, expressing the K and tt masses 

(1 + ^)Ml = 2ABo + i^ i8 

■K TT 

(1 + ^L 5 )M| = (r + l)mB + ^L s (57) 

■K TT 
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we can solve for mB and L 8 . For the latter, we find 

. U Fl 



+ 



r% 



(58) 



2 ' 8M 2 r\ - 1 

In order to evaluate L 8 , we need to know the value of r. If we use the number given in ||, 
r = 25.7 (with a 10% uncertainty), we find Lg = 1.17 10~ 3 (again fairly close to Lg(M n ) = 
1.1 10~ 3 ). However, recent work on the 0(e 2 m q ) corrections to Dashen's theorem may imply 
that this value of r should be revised downwards. Indeed, have shown that r is related to 
the strong part of the K° — K + mass difference by the following relation 

Ml 



r + 1 = (r 2 + 1) 



m d 



m,. 



M\ 



(59) 



m s -rh (Mf o - M 2 K+ ) QCD 
In this relation, || argue that R = (m s -m)/(m^-m„) can be estimated from the baryon sector 
or from uo — p mixing in a consistent way, givingfl2| R = 43.5 ± 3.2. The strong interaction 
part of the K — K + mass difference can be written as 

(M| - M 2 K+ ) QCD = M 2 K0 - M 2 K+ + M 2 + - M 2 + 6(M 2 + - M 2 ) (60) 

Dashen's theorem^] states that 5 is a quantity which is 0(e 2 m q ). Recent estimates |27[ yield 
S to be positive and in the range 0.4 < 5 < 0.8. In the apparent absence of any good reason to 
doubt the estimate of R, using this in relation (|59"D implies that the central value of r should 
be moved down to r ~ 23. This has the consequence of increasing the value of L$ up to 
L 8 = 1.43 10~ 3 . This innocent looking modification has a strong influence on the eta decay, as 
we will see. 



It is interesting to compare these results with sum-rule estimates. For L 5 [lOj gives 

F 2 
L 5 = —^ ~ 2.310~ 3 (61) 



(using Ms = M ao = 980 MeV ) which is in fairly good agreement with the above evaluation. 
A less publicized chiral sum-rule can also be derived [|n[ giving L 8 in the following form: 

L* = JL(l + t a n 2 0(l-^L)) (62) 



16M| 



where G is defined in (|15D. Using the value of Lg found above one obtains the estimate of B 
used in sec. 2. 

Let us now consider the masses of the t] and the r)'. At order four we have not only a mass 
matrix M but also a kinetic matrix K. Both matrices are real and symmetric and the matrix 
elements of the mass matrix are 
1 



M 



ii 



M 12 

M 22 



-M n 
V2. 



2r 2 + l + 2r 2 (r 2 + 2)z + 4(r 2 



ify 



(r 2 - l)M n 2 1 + (r 2 + 3) z + 2(r 2 -l)y + 2k 2 - h 



M 2 (l 



2kx) + -M v 



r 2 + 2 + ((r 2 + l) 2 + 2) z + 2(r 2 - I) 2 y + 2(r 2 + 2)(2A; 2 - k, 



(63) 
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where 

y = 4M%L 8 /Fl z = 4M 2 L 5 /F 2 (64) 

The mixing with the tt° has a negligibly small effect at this level and will be ignored. The 
matrix elements of the kinetic matrix now are 

K 11 = l + ~{2r 2 + l)z } K 12 = -^(r 2 -l)z, K 22 = 1 + ~(r 2 + 2)z (65) 

The mass squared of the 77 and 7/ mesons are given by the eigenvalues of the generalized 
eigenvalue problem: 



M n M 12 \ I V u V l2 \ I K u K l2 \ I V n V 12 \l M 2 
M 12 M 22 ) { V 21 V 22 ) " { K 12 K 22 ) { V 21 V 22 j { M4 



, t (66) 



As a result, the mixing matrix V is not unitary (instead of satisfying 4 VV=1 it satisfies 
*VKV=1). This implies that one has to introduce two mixing angles 6q and 8$ and not just 
one as is usually done (this fact was pointed out in ref . |U| ) : 

08 = !~{<f>V COS #8 + </y sin 0s) 

^8 

0o = -r(-077 sin 9 + 4>^ cos6» ) (67) 

^0 

The three parameters k%, k 2 and Mq which appear in the mass matrix can effectively be absorbed 
into two unknowns Mq (1 — 2ki) and 2k 2 — fci, which we expect to determine from fitting the 
r) and 7]' masses. It turns out that one has two different solutions, which we label as (+) and 
(— ). The requirement that the large N c expansion be meaningful allows one to eliminate one 
of the solutions. Indeed, numerically, the solution (+) corresponds to the following values for 
the parameters M and k 2 (setting k\ =0): 

M = 898 MeV k 2 = 0.12 (68) 

while the solution (— ) corresponds to 

M = 1226 MeV k 2 = -1.43 (69) 

Clearly, the solution (+) does satisfy the criterion that the contribution of the k 2 term in the 
mass matrix is roughly of the same magnitude as, say, the contribution from the L§ term. This 
does not hold for the other solution. The possibility of reproducing exactly the 77 and 77' masses 
while having only small corrections to the 0(p 2 ) lagrangian (|5T| ) was noted by |[32| . One can 
express the mixing angles in closed form 

( Mn-KnM^ Y ± -K 12 
tan6» = - 2 — '— , 6> 8 = O ± 0, tan0 = — y (70) 



KnM^-MnJ ° (K U K 22 -K 2 2/ 



1 
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as well as the extra factors 

A 8 = (K n - §^) 5 A = (K 22 - |^) I (71) 

-^22 An 

In practice, however, these expressions should be expanded linearly in terms of the 0(p 4 ) 
parameters in order to remain consistent with the 0(p A ) precision. Before exploiting this, let 
us examine the results that one would obtain at order two of the mixed expansion, i.e. starting 
from fl5"T|). The unique parameter in that case, Mq, may be adjusted so as to reproduce the r\ 



mass. The mixing angle and the r/ mass then satisfy 

tm9p = _ (2,- 2 + l)AC-3M/ & =m m^ + ^^| (72) 



(2r 2 + l)il4 2 - ?>M 2 


M 2 


= M 2 


/2r 2 + l ! V2(r 2 -1) 


V2{r 2 -l)M n 2 


rf 




\ 3 3 tan 6*p 



Numerically, this gives M„/ = 1583 MeV, which is indeed larger than the experimental mass as 



expected from the bound derived in|33|], and 9 P = —5.6°. Note that the result differs from || 
because their mass matrix includes the term k 2 which we have considered here to be of higher 
order. In the mixed expansion the mixing angle is of order 0(1) rather than 0(m q ) as in the 
chiral expansion yet, numerically, it comes out to be rather small at leading order. The ratio 
F^/F^ differs from one already at leading order and is given by F n /F v = 1/ cosOp. 

At order four, let us concentrate on the quantities which are relevant in the discussion of 
the 7/ decay amplitude that is, tan6*o and F n /F v . The former is given in closed form by eq.([70l) 
and the latter can be expressed as 

K = cos(6 -e 8 ) 

F v \ 8 cos9 K ' y ' 

As stated above, we must expand the expressions for tan# and F^/F^ linearly in terms of the 

0{p A ) parameters y, z and the mass difference M v 2 — M 2 ,. For this purpose, let us introduce 

the notations 

_, _ 3 M/ - (2r 2 + l) M n 2 

V2{r 2 - 1)M W 

One finds for the ratio F v /F v 



p = tan# P , p= ' f - i — iXa / 2 (74) 



The first two terms in the parenthesis may be recognized as those of the standard chiral expan- 
sion at 0(p 4 ) in the leading N c approximation. Here, this formula is seen to receive corrections. 
The expansion of the mixing angle, now: 

- tan# = ip(3 - ppf) + ^i(l - p 2 )(2y - |) (76) 

Let us now return to the question of the electromagnetic widths. In the effective lagrangian 
language, the new feature compared to sec. 3 is that we have a piece containing (f> 

£^2 7 = , ^ ; K U F^ { -i- + Ak, + ^^M w \ ) (77) 

F y/6{1 + 2k i) l 32vr 9 J 
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It can be seen that the amplitude for rj' —> 27 derived from (|77j) conforms to the general 
structure predicted in |34] . 



The lagrangian for <p 8 is the same as before ([24]) with the crucial difference that the t[ term 
is absent (to be more precise, it no longer contains the a! / M^ piece in (|2l|) and the rest being 
0((N C )°) does not count at this order). In order to discuss the 77 and 7/ decays we expand 0o 
and <f)g using (^). In order to be consistent with the 0(p 6 ) precision we use the mixing angles 
at 0(p A ) in the 0(p 4 ) piece of the WZ lagrangian and the mixing angle Op in the piece which 
is already 0(p 6 ). 

It is clear from ( |77| ) that we can make no definite prediction for the decay of the rf . All we 
can do is use the experimental information to constrain the combination of the parameters k\ 
and k 3 which appear in (|77|) . Once this is done, we can make a prediction for the 77 decay. We 
can recast the eta decay amplitude in a similar form as before fl27]): 

Note that only one of the mixing angles shows up in this expression. The expansions of F^/F^ 
and tan 6*o have been determined above (^|),(^|)- It is interesting to investigate the sensitivity 
of the result upon small variations of the quark mass ratio r. We will consider the two cases 
r = 25.7, as in || and r = 23 as suggested from the corrections to Dashen's theorem. In the 
first situation one obtains 

§ = 0.703, 9 = -20.4°, r(?7 -> 77) = 439 eV (r = 25.7) (79) 

F v 

This value of r is the same as the one used in the chiral expansion in sec. 3 (standard case). The 
value of the mixing angle obtained here is larger in magnitude than the one obtained in sec. 3. 
One could a priori expect that the the resulting 77 width should be larger as well but it turns 
out that this is not the case. Indeed, there is an extra ingredient in the width which is F^/F^ 
and using the mixed expansion rules we find a result which is smaller than the chiral expansion 
result F n /F v = 0.77 ± 0.03. This perhaps suggests that the uncertainty on this parameter is 
somewhat underestimated. These two differences cancel out to some extent and one ends up 
with practically the same result as before for the 77 width. Now for the smaller value of r we 
obtain 

-f = 0.710, 9 = -23.9°, T(r/ -> 77) = 529 eV (r = 23.0) (80) 

F v 

This illustrates that the uncertainty of the order of 100 eV which was obtained in sec. 3 can be 
traced, to a large extent, to the uncertainty (of the order of 10%) on the value of the quark 
mass ratio m s /m in the standard xpt- 

5. CONCLUSIONS 

In this paper we have discussed the chiral expansion of the amplitudes for n° — > 77 and 
V ~ > 77 beyond leading order. The main new contribution in this topic is the attempt made 
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to evaluate precisely, on the basis of sum-rules, the two parameters t\ and t[ which appear in 
£jg) Z . The dominating contribution, numerically, is t[ in good agreement with the fact that it 
is 0(N%) as it picks up a pole from the 7/ meson. We believe that our estimate of t[ should 
be on the same level of reliability as the sum- rule for L7. This sum-rule is, in fact, one of the 
ingredients in the evaluation. The contribution from t\ was found to amount for roughly 10% 
of the width. Its evaluation is certainly not as precise as for t[ but we argued that the sign and 
the order of magnitude should be the correct ones. An interesting peculiarity of this parameter 
is that it does not reflect the property of a particular resonance but rather of the continuum, 
which is matched to the QCD asymptotic behaviour. It would be interesting to compare these 
results with those which obtain in other approaches, in particular improved variants of the 



NJL model which seem to perform well for the standard 0(p 4 ) LE parameters [35] . On the 
practical side, we found that the 0(p 6 ) correction raises the value of the width to a level which 
is compatible with the photon-photon experimental results. The prediction has a rather large 
uncertainty, of the order of 25% which is generated by the parameter L7 and, to some extent, 
by F n /F v . We have also investigated the relevance of these two amplitudes in connection with 
the generalized %pt which embodies values of r much smaller than r 2 . It turned out that, 
even though individual factors are very different, they combine to give a result which is also 
compatible with experiment. 

One of our aims was to test the convergence of chiral perturbation theory in the anomalous 
sector. Since the lowest order result is so small compared to experiment it was by no means 
obvious how the 0(p G ) correction could manage to bring the two in agreement. This can be 
understood qualitatively if one assumes the simple rule that a 1/N C correction has roughly the 
same magnitude as a quark mass correction. The leading term in the 0(p 6 ) correction goes as 
O(N^) x 0(m q ) so, according to this rule, it is natural to expect that it could be of a similar 
size as the leading chiral order contribution which goes as 0(N C ) x 0(1). Terms of still higher 
order, on the other hand, should be much smaller. On the basis of this counting scheme it is 
possible to include the r/ into the effective lagrangian in a systematic way. We have discussed 
the question of the 77 decay from this mixed expansion point of view. In a sense, this proves 
slightly disappointing because, even in this scheme, the essential contribution to the mixing 
angle comes at next to leading order. The mixing angle turns out to be slightly larger than in 
the chiral expansion but the predictions for the width are nearly identical. There remains to 
explore whether this kind of expansion provides some non-trivial constraints in the rf sector. 
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